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Abstract 

We construct a scalar field based cosmological model, possessing a cosmological singularity characterized by a finite 
value of the cosmological radius and an infinite scalar curvature. Using the methods of the qualitative theory of 
differential equations, we give a complete description of the cosmological evolutions in the model under consideration. 
There are four classes of evolutions, two of which have finite lifetimes, while the other two undergo an infinite expansion. 
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1. Introduction 

The discovery of the phenomenon of cosmic ac- 
celeration p[] has stimulated a study of a huge vari- 
ety of cosmological models, based on perfect fluids, 
scalar fields, tachyons etc [2 . The cosmological mod- 
els based on scalar fields have a long history, being 
used for exploration of possible inflationary scenar- 
ios [3] and for description of dark energy [4]. This 
has attracted a special attention to the technique 
of reconstruction of scalar potentials reproducing a 
given cosmological evolution [5] . 

On the other hand this development of model- 
designing art has revealed cosmological evolutions 
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possessing various types of singularities, sometimes 
very different from the traditional Big Bang and Big 
Crunch. The most popular between them is, per- 
haps, the Big Rip cosmological singularity [617] aris- 
ing in superaccelerating models driven by some kind 
of phantom matter [5] . Other types of singularities 
are sudden singularities [9], Big Brake [TO], and so 
on llj. Here we would like to study the singularities 
which are close to the known Big Bang, Big Crunch 
and Big Rip singularities, but arising at finite values 
of the cosmological factor (different from zero and 
infinity as well). Similar singularities were recently 
considered in [I^ . 

In the present paper we construct potentials 
which can drive the cosmological evolution towards 
(or from) such singularities. Combining qualitative 
and numerical methods we study the set of possi- 
ble cosmological histories in the suggested models 
to show that the presence of such singularities in 
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a cosmological model under consideration depends 
essentially on initial conditions and that the same 
model can accomodate qualitatively different cos- 
mological scenarios. 

The structure of the paper is the following: in 
Sec. 2 we construct some potentials corresponding 
to evolutions with "soft" singularities. In the third 
section we analyze their dynamics. The conclusion is 
devoted to an interpretation of the obtained results. 

2. Construction of scalar field potentials 

We shall consider flat Friedmann models with the 
metric 

ds^ = dt^ - a^(t)dl^. (1) 

The Hubble parameter h{t) = a /a satisfies the 
Friedmann equation 

= £, (2) 

where e is the energy density and a convenient nor- 
malization of the Newton constant is chosen. Differ- 
entiating equation ^ and using the energy conser- 
vation equation 

e^-ih{e+p), (3) 
where p is the pressure, one comes to 



(4) 



If the matter is represented by a spatially homoge- 
neous minimally coupled scalar field, then the en- 
ergy density and the pressure are given by the for- 
mulae 



(5) 
(6) 



where V{(j)) is a scalar field potential. Combining 
equations (g]), (g]), ©we have 



and 



J,2 



3 



(7) 



(8) 



Equation ([7]) represents the potential as a function 
of time t. Integrating equation ([8|) one can find the 
scalar field as a function of time. Inverting this de- 
pendence we can obtain the time parameter as a 
function of (f) and substituting the corresponding for- 
mula into equation ([7]) one arrives to the uniquely re- 
constructed potential V{(l)). It is necessary to stress 



that this potential reproduces a given cosmological 
evolution only for some special choice of initial con- 
ditions on the scalar field and its time derivative. 

It is known that the power-law cosmological evo- 
lution is given by the Hubble parameter h{t) ~ j. 
We shall look for a "softer" version of the cosmolog- 
ical evolution given by the law 



(9) 



where 5* is a positive constant and 0<Q!< l.Att = 
a singularity is present, but it is different from the 
traditional Big Bang singularity. Indeed, integrating 
we obtain 



In 



a{t) 



S 



(10) 



a(0) 1 - a 
If t > the right-hand side of Eq. pU|) is finite and 
hence one cannot have a(0) = in the left-hand side 
of this equation, because it would imply a contra- 
diction, making In divergent. Hence a(0) > 0, 



while 



S 



a = a(0)— exp 



S 
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(11) 



This type of singularity can be called soft Bing Bang 
singularity because the cosmological radius is finite 
(and non-zero) while its time derivative, the Hubble 
variable and the scalar curvature are singular. It is 
interesting to note that when t oo both a(t) and 
d(t) tend to infinity, but they do not encounter any 
cosmological singularity because the Hubble vari- 
able and its derivatives tend to zero. 

Let us reconstruct the potential of the scalar field 
model, producing the cosmological evolution ([9]) us- 
ing the technique described above. Eq. ([8]) gives 



2 c+i 
±\l-aS t-—. 



(12) 



We shall choose the positive sign, without loosing 
generality. Integrating, we get 



m = ^3«5 ^— -, (13) 

up to an arbitrary constant. Inverting the last rela- 
tion we find 



o \ 1/2 1 

3 \ 1 — a 



2aS 



(14) 



Hence, using Eq. ^ we obtain 



aS 



(15) 
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This potential provides the cosmological evolution 
([9]) if initial conditions compatible with Eqs. (|12p 
and (|13p are chosen. Naturally, there are also other 
cosmological evolutions, generated by other initial 
conditions, which will be studied in the next section. 

3. The dynamics of the cosmological model 
with a — ^ 

In order to achieve some simplification of calcula- 
tions we shall consider a particular model, namely 
the one with the choice a = ^ . In this case 



and 



a{t) = a(0)e2^^. 



(V30/2)4 3(730/2)6- 
The Klein-Gordon equation reads 



(16) 
(17) 



+ 3<A s.gn(M^ -0^ + - 

32V354 32^/354 _ 
(V30/2)5 (V30/2)7 

This equation is equivalent to the dynamical system 



(f> = X, 



1654 
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3xsign(/i)y 2 + .^^^^^^^ 3(730/2)6 
32\/354 32V354 



(V30/2)5 (V30/2)7- 

(19) 

The qualitative analysis of dynamical systems in 
cosmology was presented in detail in [13] . 

First of all, let us notice that the system has two 
critical points: 6 = ±^,x = 0. We consider the lin- 
earized system around the point with the positive 
value of (b: 



4 o . (20) 

i = -3x sign(/i) — 5^ 32 • 35'4(^, 
v3 

where ip = (f> — 2/^/3. 

The Lyapunov indices for this system are (for h > 

0) 

Ai = -8V35^ (21) 
A2 = 4V35^ (22) 
For negative h, corresponding to the cosmological 
contraction, the signs of Ai and A2 are changed. The 



eigenvalues are real and have opposite signs, hence 
both the critical points are saddles. The universe 
being in one of these two saddle points means that 
it undergoes a de Sitter expansion or contraction, 
according to the sign of h, with the value of h — 
given by 



. 452 



(23) 



For each saddle point there are four separatrices 
which separate four classes of trajectories in the 
phase plane x, (f> corresponding to four types of cos- 
mological evolutions. 

In order to simplify the study of the dynamics let 
us note that the potential is an even function of the 
scalar field (p and that the saddle points are also 
symmetrical with respect to the x axis. Thus, it is 
sufficient to consider only one of this saddle points. 
We shall carry out our qualitative analysis taking 
into account both Fig. [U giving the form of the po- 
tential, and Fig. [21 representing the phase portrait 
in the plane 0, a;. 

First let us consider trajectories which begin at 
the moment t = 0, when the initial value of the 
scalar field is infinite, its potential is equal to zero 
and the time derivative of the scalar field is infinite 
and negative. In terms of the Fig. [1] it means that 
we consider the motion of the point beginning at 
the far right on the slope of the potential hill and 
moving towards the left (i.e. towards the top of the 
hill) with an infinite initial velocity. Such a motion 
for h > describes a universe born from the stan- 
dard Big Bang singularity. Further details of this 
evolution depend on the asymptotic ratio between 
absolute values of and at < ^ 0. If this ratio 
is smaller than some critical value then the scalar 
field does not reach the top of the hill and at some 
moment it begins to roll down back to the right. 
During this process of rolling down the scalar field 
increases, the potential is decreasing and the veloc- 
ity (j) becomes positive and increasing. However the 
universe expansion works as a friction and at some 
moment its infiuence becomes dominant causing an 
asymptotic damping to zero of the velocity. The uni- 
verse expands infinitely with h{t) — s- 0. In the phase 
portrait (Fig. [2]) such trajectories populate the re- 
gion //. This region is limited by the separatrices 
(3 and 7. The first one corresponds to the positive 
(for h > 0) eigenvalue A2, while 7 corresponds to 
Ai. If the ratio (j)/4> has the critical value, then the 
scalar field reaches asymptotically the top of the hill 
of the potential, meaning that the universe becomes 
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asymptotically de Sitter: in the phase portrait it is 
nothing but the curve 7. 

When the ratio introduced above is larger than 
the critical one we encounter a different regime. In 
this case the scalar field passes with non- vanishing 
velocity the top of the hill and begins to roll down 
in the abyss on the left. The absolute value of the 
negative velocity </> is growing, while the potential 
becomes negative and at some moment the total en- 
ergy density of the scalar field vanishes together with 
the Hubble parameter h: this means that the uni- 
verse starts contracting. This contraction provides 
the growing of the absolute value of the velocity of 
the scalar field and the kinetic term again becomes 
larger then the potential one. Moreover, both the 
terms in the Klein-Gordon equation increase the ab- 
solute value oi (p < 0. One can easily show that 
the regime in which the time derivative (f) becomes 
equal to —00 at some finite value of (j) is impossi- 
ble, because it implies a contradiction between the 
asymptotic behaviour of different terms in Eq. ()18p . 
Thus, the universe tends to the singularity squeez- 
ing to the state with the value of equal to zero and 
an infinite time derivative (p. To understand which 
kind of singularity the universe encounters, we need 
some detail about the behaviour of the scalar field. 
Let us suppose that, approaching the singularity at 
some moment to, the scalar field behaves as 



(24) 



where < /i < 1. Then the first and second time 
derivatives are 

Pit) - -fiMto-tr-\ m = Ai(Ai-i)0o(io-i)^"'. 

(25) 

The potential behaves as 



V 



20485"* 



(26) 



To have the Hubble variable well defined we re- 
quire that the kinetic term is larger than the abso- 
lute value of the negative potential term (|^^ . i.e. 
2fi — 2 < — 6 fi, or fi < J. Now two opposite cases 
may hold: (i) the friction term in the Klein-Gordon 
equation could dominate the potential term or (ii) 
the opposite situation. For (i) to be the right case, 
one should require 2/i— 2 < — 7/ior/i < |.Inthissit- 
uation the asymptotic behaviour of the second time 
derivative of (p should be equal to that of the fric- 
tion term, or, in other words /i — 2 = 2/i — 2, that is 
/i = 0, which obviously is not relevant. Thus we have 
to consider the range | < /i < i. In this case the 



potential term should be equal to the second time 
derivative of (p, which implies: 



and 



yo = 4a / — . 



(27) 



(28) 



Substituting the values of fj, and (po into the expres- 
sion for h{t), we obtain 



h{t) = - 



S 



(29) 



Thus, we see that the singularity we are approaching 
is of the soft Big Crunch type. 

In the phase portrait (Fig. ^ these trajectories 
occupy the region /// limited by the separatrices 
7 and 5. The cosmological evolutions run from the 
Big Bang singularity to the soft Big Crunch one. 
However, the Fig. [5] is not sufficient to describe the 
complete behaviour of the universe under consider- 
ation, because at some moment the Hubble variable 
changes sign and we should turn to the Fig.[3l giving 
the phase portrait for the contracting universe h < 
0. Note that increasing the velocity ip with which 
the scalar field overcomes the top of the hill, implies 
delaying the time when the point of maximal expan- 
sion of the universe is reached. 

The third regime begins from the soft Big Bang 
singularity, when the scalar field is equal to zero and 
its time derivative is infinite and positive. In Fig. [1] 
that situation is represented by the point climbing 
from the abyss to the top of the hill. If the velocity 
term is not high enough, at some moment the field 
stops climbing and rolls back down. During this fall 
the Hubble variable changes sign and the universe 
ends its evolution in the soft Big Crunch singularity. 
The corresponding trajectories belong to the region 
IV of our phase plane, bounded by the separatrices 
S and a. The universe has its finite life time between 
the soft Big Bang and the soft Big Crunch singu- 
larities. The situation when the scalar field arrives 
exactly to the top of the hill and stops corresponds 
to the separatrix a. 

The fourth set of cosmological trajectories is gen- 
erated by the scalar field climbing from the abyss 
and overcoming the top of the hill with the subse- 
quent infinite expansion: the scalar field is infinitely 
growing and the Hubble parameter tends to zero. 
These trajectories occupy the region / and our origi- 
nal cosmological evolution ([9]) belongs to this family. 
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4. Conclusion 

Let us sum up our results. Wishing to describe 
a cosmological evolution beginning from (or end- 
ing in) a singularity characterized by a nonvanish- 
ing initial (or final) cosmological radius and an infi- 
nite value of the scalar curvature due to the infinite 
value of the Hubble parameter, we have constructed 
a scalar field potential, providing such an evolution. 
Then, using the methods of qualitative analysis of 
the differential equations, we have shown that the 
proposed model accomodates four different classes 
of cosmological evolutions, depending on initial con- 
ditions. Numerical simulations have confirmed our 
predictions. 

The main results of this work are: (i) the realiza- 
tion of a concrete cosmological model with a scalar 
field, where finite cosmological radius singularities 
are present and (ii) the complete description of all 
the possible evolutions of this model depending on 
the initial conditions. It is important to remark that, 
given the fixed scalar field potential, one has differ- 
ent types of evolutions encountering different kinds 
of singularities. The trajectories, belonging to the 
regions / and II have an infinite time of expansion, 
while the trajectories belonging to the regions III 
and IV begin and end their histories in the singu- 
larities and have a finite lifetime. 

Let us note, that there are some studies |9|ll|12j 
devoted to the general analysis of various kinds of 
new cosmological singularities. In this work we were 
not looking for an exhaustive classification of differ- 
ent possible cosmological models possessing singu- 
larities, but rather we wanted to study in a complete 
way a particular cosmological model, having some 
interesting properties. 
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Fig. 1. Plot of the potential V((/<) as given by Eq. JTtJ. Here 
only the positive t/i-axis is shown, which is enough to under- 
stand the behaviour thanks to the parity of the function. We 
have also delineated the four different dynamical behaviours 
of the scalar field, clearer to understand taking into account 
the phase portrait (see Fig[2]|. 





Fig. 3. Phase portrait for the dynamical system under discus- 
sion (Eq. I I19I I) with h < 0, i.e. describing evolutions of the 
universe characterized by contraction. As for the h-positive 
case, we can see the four separatrices of the saddle point 
and the corresponding four regions of different dynamical 
behaviours of the trajectories. 



Fig. 2. Positive </i-axis phase portrait for the dynamical sys- 
tem II19I I with h > 0. The four separatrices of the saddle 
point (a, /3, 7, (5) individuate four regions {I, II, III, IV) with 
different behaviours of the trajectories, as explained in the 
text. 
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